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ABSTRACT 

Based upon the bubble nucieation and site deactivation 
mechanisms proposed by Kant, a dynamic model was developed to 
study the stability of nucieation sites in the heterogeneous 
pool boiling of liquids on a singie in a natural heated surface. 
This model predicted the waiting time of a bubble for given 
boiling conditions in agreement with the available experimental 
data. The results showed that the diameter and the depth of the 
nucieation site actually determined the stability and other 
parameters in pool boiling were not important. 

An equation was developed for determining the thickness of 
the thermal layer on a heated surface. Computed values of the 
thickness were found to be in agreement with the values availa- 
ble in literature for water and pentane. 

Static single bubble nucieation model of Kant was extended 
to incorporate a generalised temperature distribution for the 
thermal layer, proposed here. This model was used to predict 
the nucieation superheats and the size range of active sites on 
a surface. The results obtained were better than those of 


other authors 
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CHAPTER 1 


INTRODUCTION 


1.1 POOL BOILING 

In a normal pool boiling situation, the surface of the 
container is usually at a temperature higher than that of the 
liquid contained therein and this liquid may be at or near its 
saturation temperature. The heat transfer takes place from 
the surface to the liquid mostly through a process of phase 
change of some of the liquid into the vapor state. The 
applications of such boiling heat transfer are many, not only 
in apparatuses whose primary purpose is to vaporize a liquid, 
such as boilers, but also in situations where it is desired 
to remove heat from a surface at a high rate with the lowest 
possible surface temperature, such as, nuclear reactors. 

There are several regimes of pool boiling, like nucleate 
boiling, transition boiling and film boiling. The most 
desirable and the important regime is that of nucleate boiling 
in' which bubbles are formed wherever there are nucleation 
sites such as tiny pits or scratches in the surface. The 
bubbles transport the latent heat of the phase change and also 
increase the convective heat transfer by agitating the liquid 
near the heating surface* Nucleate boiling is essentially a 
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cyclic process which consists of bubble formation, bubble 
growth, and its departure. The nucleation sites have to be 
active or stable for the process to continue. The subject of 
this thesis is to model the bubble nucleation phenomena to 
predict the stability of nucleation sites. 

Some of the terms used in the following chapters are 
defined here. Heterogeneous pool boiling refers to the boiling 
occurring on a surface or on a single cavity in the surface 
submerged in a pool of liquid. If the liquid is at the 
boiling point, the term saturated pool boiling is used and if 
the liquid temperature is lower than its boiling point, the 
term sub-cooled pool boiling is used. There is a small layer# 
of liquid adjacent to the boiling surface, called the thermal 
layer, in which most of the temperature drop occurs from the 
surface temperature to the bulk liquid temperature. 

1.2 OBJECTIVES 

1. Develop a general equation for the thermal layer thick- 
ness . 

2. Propose a generalized temperature distribution within the 
thermal layer and predict nucleation superheats. 

Develop an expression to predict the size range of active 
site on a surface. 


3 . 
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4. Develop a dynamic single bubble nucleation model for the 

heterogeneous pool boiling (KP3) of a liquid to predict 
the minimum wall superheats for the stability of 
nucleation sites in a surface. 

1.3 ORGANISATION OF THE THESIS 

Chapter 2 gives the general background of the subject and 
summarizes all the equations which are used later in the text. 
Chapter 3 gives the expression for the thermal layer thickness 
and those for the generalized temperature distribution within 
the thermal layer and the size range of active nucleation sites. 
The results obtained are compared with the predictions of other 
workers. Chapter 4 gives the dynamic single bubble nucleation 
model for HPB. The predictions of the model are discussed and 
compared with the predictions and data of other authors. 

Chapter 5 gives the conclusions followed by references. 



CHAPTER 2 


GENERAL BACKGROUND 


2.1 INTRODUCTION 

This chapter deals with the topics that are referred to 
later on in the thesis. A comprehensive review in Boiling Heat 
Transfer is not attempted here. 

2.2 FUNDAMENTAL EQUATIONS 
1. Laplace Equation 

For a curved interface, the pressures in the liquid, p^ , 
and that in the vapor, p , are related by the Laplace equation 



where R^ and P^ are the principal radii of curvature of the 
interface and a is the surface tension of the liquid. For a 
spherical interface = Pq 



( 2 . 2 ) 


2. Clausius-Clapeyron Equation 

V. 

The pressure difference (p v ~p^ ) can be related to the 

equilibrium vapor superheat (T **T ) through this equation and 

G S 
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this relation is given by 


dp 

dl 



(2.3) 


where P v and represent the density of vapor and enthalpy 

of vaporization respectively. 

Equation (2.3) upon integration along the p-v saturation 
P hr 

curve assuming — = constant and T = T s (Rohsenow, 1970), 
yields 


201 

T — T = — J 

e s R. P h. 

x v fg 


(2.4) 


2.3 THERMAL LAYER THICKNESS 

The thermal layer thickness, & , (defined in Section 1.1) 
plays a very significant role in the case of nucleate boiling. 
It also plays a very important role in determining the size 
range of effective nucleation sites. There are some empirical 
relations available in the literature to estimate the value of 
6 + . Most of these relations, relate 6 + to just one or two 
properties of the liquid like k<j » v , x , a and the heat trans- 
fer coefficient, h , ana are based upon limited and poorly con- 
ducted experiments. 
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Lippert and Dougall (1968) for the pool boiling of water, 

Freon-113 and Methyl Alcohol showed that the correlations for 

5 £ could be divided into two distinct regimes. For low heat 

-2 —1 

fluxes (h <_ 3970 Vv m K ), for ail the three liquids 

5 t = C k /h (2.5) 

where is the thermal conductivity of the boiling liquid 
and the constant C was taken as 1.65. 

For higher heat fluxes the correlation was 

6 t = (C/h) l/2 * 6 (2.6) 

where the value of C depended on the particular liquid. 

According to Marcus and Dropkin (1965), for low heat 
fluxes (h < 3970 V/ K 

6 1 = 1.57 k./h (2.7) 

and for high heat fluxes (h > 3970 'fl K 2 ) 

6 1 = 0.0177 h“ 0 * 5 metre (2.8) 

Han and Griffith (1965) expressed 6 ^ in terms of the 
waiting period, t as 

6. = (nat ) ljf2 
t v w 


(2.9) 
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where a is the thermal diffusivity of the boiling liquid, and 
t is the waiting time of the bubble (defined in Section 2.5). 

For very high wall superheat, where the waiting period 
is very short, the thickness, was given by 

«t = [* “(t w + t ub )] 1/2 (2.10) 


where t ^ is the unbinding period. 

As the waiting time increases the thermal layer cannot 

increase indefinitely. It will be washed off by the natural 

convection currents of the fluid as it grows beyond the 

thickness of the natural convection layer, 6 • This means 

nc 

that the upper bound on o^_ can be written as 


max 


= 5 


nc 


( 2 . 11 ) 


The lower bound on 6 ^ can be determined by using an h value 
for nucleate boiling from iiohsenows* (1952) correlation 


(T -T > 
L_w s ' 

‘ h 

fg 


= r • 


sf 


q 

ITT 


f g 


a 0.33 

{ * — - } j 

g( p-, -■ p v ) 


c 

(- 


1 p. m+1 

( 2 . 12 ) 


where q , , p and represent the heat flux, specific 

heat, viscosity and a constant for the surf ace-liquid combina- 


tion . 



8 


Leont'ev and Kirdyashkin (1966) gave the following expre- 
ssion for predicting 6^, 


2 *75(- 


\ l/3 


(2.13) 


where a, v and A. represent the thermal diffusivity , kinema- 
tic viscosity and coefficient of volume expansion. 

2.4 TSiiPERATURE DISTRIBUTIONS IN THE THERMAL LAYER 

Following temperature distributions (TD) have been 
suggested for the thermal layer. 

1. Linear Temperature Distribution (LTD) 

The LTD is given by 




(2.14) 


2. Nonlinear Temperature Distribution (NLTD) 

Leont’ev and Kirdyashkin (1966) approximated the TD 
by the following equation 


£ = 1 - 2r| 4 - 2n 2 3 - ti 4 


where 



V 



(2.15a) 


and 


(2.15b) 
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3. Marcus and Drop kin Temperature Distribution (MDTD) 

Marcus and Dropkin (1965) suggested the TD as % 

T - T 

= l - --X 0 < y < a 6. (2.16a) 

w b 0 1 

T - T 

^ = d ( 7 Xj-n a6 < y < 6 (2.16b) 

w b °t 

where a,d and n are constants. 

2.5 STATIC BUBBLE NUCLEATION MODELS IN HPB 

The presence of thermal layer on a surface makes the predi- 
ctions of incipient wall superheats difficult, compared to a 
situation where the boiling liquid is uniformly superheated. The 
latter situation corresponds to = «>. 

As mentioned earlier, the formation of a bubble is a cyclic 
process , the time period between the departures of two succe- 
ssive bubbles is the total period of the bubble cycle. This 
period is comprised of two periods, namely, the waiting period, 
t w , and the growth period, t . As the names imply, the former 
refers to that period during which conditions of heating are 
reset to produce a visible vapor nucleus on a cavity in the sur- 
face following the bubble departure. The latter, i.e., t , 

g 

refers to the period during which the above mentioned visible 
nucleus grows to its maximum departure size. The reciprocal of 
(t w + t ) gives the frequency of bubble formation. Several 
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definitions of these terms appear in the literature (Hsu and 
Graham 5 1961; Han and Griffith, 1965; Shoukri and Judd, 1978). 

Many static models have appeared in the literature which 
basically have attempted to develop expressions for the wall 
superheats required for the bubble nucleation on a surface. 

The prominent ones are briefly described below. 


Most of these assume a LTD 'within the thermal layer. For 
known wall superheats all these models can give a size or the 
size range of active nucleation sites. 

1. Griffith and Wallis( 1960) observed that equation (2.4) 
predicted the required superheat only for a uniformly super- 
heated situation and this value was much lower than the actual 

» 

v/ all superheat required for nucleating a bubble on a surface in 
HPB. Using artificial conical cavities of 51 pm diameter, 
under HPB, they observed that nucleation commenced at a wall 
superheat of 11K rather than the value of 1.6 K predicted by 
equation (2.4). 


2. Hsu (1962) proposed the following equation to predict the 
size range of active sites 


( R 1 

v c ; max * min 


TO 7” 


[1 + £ 1 - 


8C b g T jT w -T b ) 

6 t P v h fg< T w- T s> 


r 1 1/2 1 


,, • (2.17) 

“ 

where C and CM are constants depending upon the bubble-surf acft 
a D ~ ^ 

contact angle and cavity geometry. The above relation gives 
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the nucleation superheat; A , for saturated pool boiling 

(T b = T s ) as 

AT n = ( V T s ! = 2 V^- { < 2c aV 6 P - 1 > 3 (2.18) 

From equation (2.17) one can deduce a stability criterion for 
sustained boiling to exist at a given nucleation site and it is 
given by 


T 

w . 
nun 

where 


T c + A h + [A h + 2 (T c -T h ) } ] 


1/2 


^ P V * fg 5 t 


s A b ; 


(2.19) 


( 2 . 20 ) 


3. Han and Griffith (1965) gave the following expressions for 
the size range of active sites for a hemispherical nucleus 
sitting at the cavity mouth 


R 


5 t ( V T s> 


'max * mm 


3 ^vv 


[i + ci- 


12aT s< T w- T b) , 

P v h fq 6 t<V T S )2 


1/2 


( 2 . 21 ) 


For T, 


V 


12 a T 


AT n = 


3R. 

P v h fg S t f 1 - (“a" - 1 ) 2 ! 


( 2 . 22 ) 


which can also be considered as the minimum wall superheat 
required to emit bubbles from a cavity of a given size. 



12 


4. Howell and Siegel (1967) developed a model for a hemi- 
spherical nucleus sitting at the cavity mouth. They used the 
concept of evaporation and condensation at the 1-v interface. 

The superheats due to them are given by 

R c < 6 1 (2.23) 

R c > 6 t (2.24) 

where (T - T ) is given by equation (2.4) with R- = R . 

5. Kant (1983) extended the nucleation model due to Howell and 
Siegel (1967) , but neither assumed a hemispherical vapor nucleus 
at the cavity mouth nor a LTD. His model equations which were 
developed for the three TD in the thermal layer namely the LTD, 
NLTD and MDTD are given below. 

LTD : 

3 (T e -T b ) 

(L-TJ > ” — — H < 5 1 (2.25) 

- R_ • R t 

[3(1 - ™ cot p) - - c ~~ ] 

6 t 2 5 t 

2BR 

( VV2T t -‘VV 


AT 


iVlJs ) 

n _ - R — 

(1 - -JL. ) 
2 6 , 


2R 

AT n = < T e - V 


H > 5 t (2.26) 
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NLTD 


‘W ->■ 


[i 


{ T - T ) 

v x e b ; 


1 , 3R Cx3 

2 ’'“ST 5 

u 


A (?^) 4 ] 


H < 6 t (2.27) 


10. 


~> 3 BR c < W 


H > 5. 


(2.23) 


JVi 


DTD 


( T w- T b> 


BR (T - T. ) 
c e b 


2 


6 t [(a 


a 


£ ql-n-. 

l-n / tl 6,° J 


2”) + (tl") ^ t* 


1-n , *i 

a } J 


H < 6 


(2.29) 


<W a 


BR ( T 
c e 


V 


« t [(a 


2 , , 

a n & /. l-n\-| 

5“ ) t X^T ^ i " a '- 1 


H >6. 


(2.30) 


where B 


_ I t- cos 3 


sm b 


(2.31) 


and p = angle which the bubble surface makes with the boiling 
surface (shape angle). 

The predicted wall superheats for the NLTD were higher 
than those for the LTD ana MDTD. Besides the effect of TD, 
his model also incorporated the effect of the shape angle, p. 
Three values of p namely £3 = n/2, (te/ 2 + ©) and G were used. 
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whore 0 was the dynamic contact angle for the surface-liquid 
combination. The value of p = n/2 corresponded to the hemis- 
pherical shape of the bubble, but not the other two values. 

The values = © and (tc/2 + 9) were shown to be important with 
reference to his dynamic bubble nucleation model (discussed 
in Section 2.6). 

2.6 DYNAMIC BUBBLE NUCLEATION MODELS IN HPB 

All static models of nucleation superheat assumed the 
presence of a vapor nucleus on the cavity mouth and the 
nucleation criteri-on was satisfied when the 1-v interface 
was exposed to the bulk liquid. Such models say nothing 
about the mechanism of nucleation and cite deactivation. The 
vapor nucleus grows when the thermal layer is re-established. 
The site is considered active if the waiting time is finite. 
Many dynamic models have been proposed to simulate the bubble 
nucleation and site deactivation mechanisms, which are briefly 
described below. 

1. Bankoff (1959) modeled the bubble nucleation problem as 
viscosity controlled cavity penetration due to vapor 
condensation. He did not predict the reversal of the 
motion of 1-v interface. 

2. Marto and Rohsenow (1966) neglected the interface motion 
in their energy balance and derived a dimensionally 
incorrect equation to predict the interface motion. 



Fabic (1964) used both advancing and receding contact 
angles? the data for which are not readily available. 

The superheat predicted by him is likely to be much 
lower than the actual wall superheat needed for 

nucleation . 

Singh et a i. (1975) developed a model for the stability 
of a nucleation site considering transient heat flux, 
inertial a nd viscous effects of the liquid. The time 
required for the interface penetration was taken to be 
zero which ± s n ot in conformity with the experiments. 

Kant (1983) based upon his boiling nucleation experi- 
ments on single sites in the pool boiling of organic 
liquids and their mixtures, threw light on the Bubble 
Nucleatio n anc j site Deactivation Mechanisms. Based upon 
his observations, he proposed a dynamic single bubble 
nucleation model which predicted the waiting time of a 
bubble for given boiling conditions. This model took 
into account the flow of liquid into the nucleation site 
and the time variation of the surface temperature through 
an additional model proposed for the same. Other feature 
of his model which were not considered in the earlier 
models for predicting t , were the condensation and eva- 
poration a t the 1-v interface, pressure losses during the 
flow and the modeling of flow within the cavity. Both 
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Hagen-Poiseuille and Suddenly Accelerated Flows were 
considered within the cavity. The predictions of the model 
were in - agreement with the experimental observations and the 
data. 

The above model established an implicit site stability 
criterion for a given cavity diameter, degree of superheat 
and degree of subcooling and a given surface-liquid combina- 
tion. His other criterion was directly based upon his obser- 
ved, site deactivation mechanism and was explicitly expressed 
as a limitation on the vapor volume within the site which 
predicted a particular value of (L/D c ) ratio for the site to 
be stable. The computation of this ratio made use of the 
model described above. The predicted ratio was in agreement 
with that for cavities in a natural surface. For initial 
penetration of the liquid, a minimum aspect ratio (L/D ) 
purely as a function of © was suggested. 



CHAPTER 3 


THERMAL LAYER THICKNESS AND NUCLSATION SUPERHEATS IN POOL 

BOILING 

3.1 INTRODUCTION 

The thermal layer thickness, 5^, as defined in Section 1.1 
has been used extensively in the pool boiling literature, but 
no proper method exists for its determination. Most of the 
static models for bubble nucleation use this parameter in their 
expressions, but its value ultimately is computed from some 
empirical expressions suggested for water and a few organic 
liquids based upon limited and rather unreliable experimental 
data. These expressions model 6^ as a function of liquid 
thermal conductivity and heat transfer coefficient. Many authors 
have also talked about the upper and lower bounds of £ ^ corres- 
ponding to the free convection and pool boiling situation res- 
pectively. An attempt is made here to develop a general expre- 
ssion for 6.J. and based upon its value predict the boiling nuclea- 
tion superheats and the size range of active sites. 

3.2 MODELLING OF THERMAL LAYER THICKNESS 

Nucleate boiling on a metallic surface is dependent upon 
the fluid properties, heat flux, degrees of superheat and sub- 
cooling and the surface roughness. It is postulated here that 
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the existing expressions for 6., in the literature, are too 
simple to be realistic for any given boiling conditions and all 
the aforementioned parameters should be included in developing 
an expression for it. This is done here using the 'Dimensional 
Analysis'. The degree of subcooling is not considered here. 


As mentioned above, <5^ can be written as a function of 
fluid and vapor properties, heat flux, <~ , degree of superheat, 
AT, cavity diameter, D , and surface roughness, S n . Then, 

6 t = fn > T w - T s > P V ,P 1 » h fg * a » D C 5 » k 1 » C 1 » S b) 

( 3 . 1 ) 

The dimensions of the different quantities are given in 
Appendix B . 

Using the Buckingham's n theorem and taking AT,P. , D , q 

1 c 

as the repeating variables, one can obtain the following eight 
it's or the dimensionless groups. 


rc, = 


5 t /D c 


= 


TCo = 




= 


V P 1 

2 , * 

P h- D / a 
v fg c / 4 

0 D r / q 


it , 


= 


^ / ( P 1 q ) 1 ^ 2 
k 1 p ^ 2 D 2 c AT/ q 3/2 
P 1 c 1 AT/ q 


( 3 . 2 ) 


*8 ~ S f/ D C 



19 


Using equation (3.2), equation (3.1) can, therefore, be written 
as 


(3.3) 

(3.4.) 

~ % 4 % Q ~ a s r/q (3.5) 

and 

ft 2*2 

n 4 = 7t 5 n 6 = k^ |i 1 D c Al/q (3.6) 

3 

Using equations (3.4) through (3.6), equation (3.1) can be 
written as 


it f n (ti^j ^ 3* ^ 4 5 ^6 J ^7 * ^ g ) 


-* 


Defining '^3 and % A as follows : 


J akob 


it 


2 ^3^7 “ ^ P 1 °1 AT// P v h fg-^ ” J a ^Number^ 


■>t -'f 


* , 


= fn(it“, ti 3 , it 2 » ^4) 


oS- 

fn( J , -ri 
a q 


p v h D c AT , 


’ p i ’ 


• 2 

q 


D 


„ h ,°% e 2 P v®o . , , . > 

{ ~T ] ( p 7 ) 1 ’ ) 


3 - 6 p 1 D c AT ' S 4 
.2 " 
q 


(3.7) 


where C s ^. - constant depending upon surface-liquid combination. 


Taking logarithms on both sides, yields 

V 

v 


e 1 log(J a ) + e 2 log(-^-B) + e 3 log(-~) + s 4 log( 


k 1 “1 D c AT . 

:7 ) 

q 


D c C sf 


= log 


(3.8) 
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The indices e^, and can be obtained by solving 

equation (3.3) if atloast four experimental values of 6 ^ and 
corresponding values of other parameters are available for a 
boiling liquid. 

3.3 BUBBLE NUCLHATIOH MODEL FOR POOL BOILING 

Kant's (1983) 'Static Single Bubble Nucleation Model' 
for determining the wall superheat required for the incipience 
of bubble formation is extended here to incorporate a generali- 
zed temperature distribution, within the thermal layer, proposed 
here. The model takes into account the condensation and evapo- 
ration at the 1— v interface. If y represents the distance mea- 
sured from the surface into the body of the liquid, then at 
y = 0, T = T w and at y = 6^, T = T, 0 j y = 5 Q denotes the posi- 
tion where T = T ♦ Evaporation at the 1-v interface will occur 
in the region 0 <_ H <_ 6 e and condensation in the regions 
< H < 6 1 and H > 6 ^ . 

The following three cases are considered (Kant, 1983) : 

Case (i); 0<H<_6 e 

Fig. 3.1 shows a vapor nucleus of height, H, on a surface 
at temperature T . Using polar coordinates (r, 0 ,¥ ), and 
assuming axial symmetry, the heat transfer to the bubble nucleus 
of radius by evaporation can be written as 




Fig. 3.1 Vapor nucleus at a cavity of radius R c . 
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O 

= 27th R£ / (T-T e ) sin 0 d0 

P 


(3.9) 


where 

p = the angle which the bubble surface makes with the 
boiling surface 

T = liquid temperature within the thermal layer 

T = equilibrium vapor temperature 
0 


Case (ii) : 5 < H < 6a. 

0 — u 

Considering both evaporation and condensation at the 1-v 
interface, the net heat transfer over the bubble surface, in 
this case as well, is given by equation (3.9). 

Case (iii) : H > 6^. 


Again taking into consideration the evaporation and conden- 
sation at the 1-v interface, the net heat transfer over the 
bubble surface is ,given by 


Q 2 = 27th [ / (T~T e ) sin 0 d0 - ( T e "T b ) (1+cos )] (3.10) 

P 

In both equations (3.9) and (3.10) R, can be replaced, for a 
spherical bubble, by 


K b = R c ,,/sin ^ 


Height, H, of the bubble nucleus is given by 



1 + cos 3 \ 
'“sin p ' 


(3.11) 


(3.12) 
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From Fig. 3.. 1, 0 ^ corresponding to y = 6 ^ can be expressed 

as 

cos 0. = cos p - ( 6 + sin p/R ) (3.13) 

o *u c 


For the three cases mentioned above, the net heat transfer 
rates and Q 2 can be determined if the temperature distri- 
bution within the thermal layer is known. 

Let the temperature distribution within the thermal layer 
be expressed as a polynomial 


t ,,234 

£ = a Q + a 1 r) + a 2 r) + a 3 r) + a 4 r) 

where 

T ~ T b y 

£ = =r and T) = --- 

T w“ T b 6 t 

From Fig. 3.1 


(3.14a) 


(3.14b) 


y = R b (cos p - cos 0) 


(3.15) 


Substituting equation (3.14) into equations (3.9) and (3.10), 
using equation (3.15) and replacing (l+cos p)/sin p by 3 yields 


2'chRt 


a, R 


x B x (T , -Tj [a n + o' 


Q 1 “ sin p " * ^w~ x b' L °o T 2 § 


V'c u 


+ 


t 


n 2_,2 

a 2 d c° 

3 67 


a Q R^3 3 a 4 R^3 4 . 2ithR?B( T^-Tj 


4 * 


_3“c 
4 6 : 


+ 


5 6 


4 


'] 


c 

sin 3 


e A b' 


(3.16) 
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2rchR n a, 

Qs - — j£ [( V T b ) a t ^ £a 0 + ^ 
sin p c 


a 2 *3 

+ + yr~ 

o 4 


*4 


- (T e -T b ) (1+cos p)] 


(3.17) 


For the nucleus to grow both > 0 and > 0 yields (see 
equations (3.16) and (3.17)) 


( T ’*'“ T b) 2 


W 


< V T b) 

'”™ - a— a- - 23 a *4 

a,R B a 0 R“B a^R^B auVB 
r i i c . z c o c . c “i 

ta o + rrr 772™ + ~t“ + 

t 3 5 4 6^ 06^ 


(3.18) 


l V T h> 


R c B(T e - T b> 


S-tK + 7 


4” 


4 -] 


(3.19) 


Equations for a hemispherical nucleus can be obtained from 
equations (3.13) and (3.19) by substituting B = 1 corresponding 
to p = 7t/2 . 

3.4 SIZE RANGE OF ACTIVE SITES 

The equations developed' in the preceding section yield 
wall superheat required for bubble nucleation as a function of 
(T -T b ) , R q and 6 + . For a given R c , (T e ~T b ) can be calculated 

”T 

from equation (2.4) and $ ^ can be computed from equation (3.3). 
If nucleate boiling superheat is known for a surface, then the 
size range of active sites in it can be determined from equa- 
tions (3.18) or (3.19). 



25 


3.5 PREDICTIONS 

Based upon the equations developed in the preceding 
sections, computations were made for the thermal layer 
thickness, nucleation superheats and the size range of active 
sites which are discussed below. 

3.5.1 Thermal Layer Thickness ( 6 ^) 

In order to determine 6 ^, the indices in equation (3.7) 
were determined using the available values of 5 ^ for pentane. 
These indices are given in Table 3.1. 

Table 3.1 
£ ' s for Pentane 


s 


1 


£ 


2 


e 


3 


£ 


4 


1.5 


1.9 


-1.3 -0.5 


Using these indices, 6 was computed for Pentane from 
equation (3.7) for given values of AT and q, which were found 
to be in good agreement with the reported values as shown in 
Figures 3.2 and 3.3 for pentane and water respectively. The 
values of 6 ^. due to Lippert and Dougall (1968) were computed 
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Fig. 3.2 Variation of thermal layer thickness with heat transfer coefficient for 
pentane. 
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from 6^ = 1.65 loj/h, whereas those due to Marcus and Dropkin 
(1965) were computed from 6^ = 1.57 /h for h <i 3970 W m^K -1 
and <5^ = 0.0177 h for h > 3970 W m - ^K~^. Hsu (1962) used 
the data of Clark, Westwater and Strenge (1959) and modeled it 

/h. The values of 6^ computed from equation 
(3.7) are also shown in the figure. It can be seen from it 
that the predicted values of 6 ^ are in good agreement with the 
values of Hsu (1962) with the values of Hsu (1962), which led 
to the conclusion that the simpler modeling 5^ = k^/h was just 
as good as the modeling presented by equation (3.7). The 
parameters q, AT and k^ seem to have more pronounced effect on 
6-j. than any other variable. But as the available data on 6 + 
is rather unreliable, the author strongly feels that this com- 
parison is inconclusive. The real comparison can be made only 
when realistic data on 6^ become available. It may be noted 
that if the values of Lippert and Dougall (1968) and those of 
Marcus and Dropkin (1965) are divided respectively by 1.65 and 
1.57, the resulting values will be found in good agreement with 
the values of Hsu (1962) and those presented here except the 
first value for Pentane and the last for water which were 
calculated using 6 t = 0.0177 h" 0,5 . 

3.5.2 Temperature Distribution in the Thermal Layer 

Three TD' s namely LTD, NLTD and MDTD for the thermal layer 
were discussed earlier. Only for the NLTD, predicted values of 


simply as 5 t = k| 
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AT^ were higher than those for othor TD's (Kant, 1983). 

The reported experimental values for bubble nucleation 
on a surface are still higher than the predicted values. Hence 
it was attempted to propose e TD in the thermal layer 'which 
would bring a closer agreement with the experimental data. The 
proposed TD is given by equation (3.14a) and the constants are 
given in Table 3.2. 

Table 3.2 

Constants in the Temperature Distribution 


a o 

a l 

a 2 


a 3 

a 4 

1 .0 

-2.8 

1.0 


2.0 

-1.0 


3.5.3 Bubble Nucleation Superheats ( AT n ) 

The model presented in Section 3.3 related AT to D , 

n c 

6,, (T -T ) and the constants of the T D in the thermal 

"C 6 S 

layer. For a given liquid boiling on a surface having a 
cavity of diameter, D c , AT n could be computed if 6^ was known. 
In the results presented below, 

6^ was computed for different values ofA(i.e., different 

values of q and (T ~T_) ) . These values of 6. were then used 

W S X 
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to compute At^ for given D c and shape angle, p. 

The effect of increasing 6 t , keeping other factors 
unchanged, was a decrease in the value of AT . As mentioned 
earlier, three values of j3 , namely, ;3 = 0, 71/2 and tc/ 2 + 9 
were used to calculate At^. These are plotted in Fig. 3.4 
showing the variation of A T n with 6 t f or aforementioned 
values of p. It can be clearly seen from these curves that 
predicted Aj n are higher from p = 0 than those for p = tu/ 2 
and ti/2 + 0. The author strongly feels that as p = 0 is a 
characteristic of a surface liquid combination while other 
values are not, the former should be considered for calculating 
AT n - Shape angle, p, as defined earlier will keep on changing 
continuously right from the moment a 1 -v interface can be seen 
on the heated surface to the moment of bubble departure. For 
these reasons, the assumption of other authors of a hemispheri- 
cal vapor nucleus sitting on a cavity mouth, is unreasonable. 

The values of AT^ were also computed for various values of 
D c> keeping unchanged. The constants a ± (i =0,1, 4) 0 f 

the TD given in table were used. The trend of the predictions 
as shown in Fig. 3.5 is similar to those of Howell and Siegel 
(1967) and Kant (1983), but the predicted values of At are 
somewhat higher than thope of other authors. 



Saturated Pool Boiling of Water 
on a copper surface 

Dynamic Contact Angle 0= 22 fl 



Thermal Layer Thickness, S\ , jjm 

Fig. 3.4 Variation of nucleation superheat with thermal layer thickness 
for water. 
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3.5.4 Size Range of Active Sites 

The equations (3.18) and (3.19) in the preceding section 
related (T W "’T S ) to (Tg-T^), R c and 0. For a given boiling 
liquid/ R c and 8, it is easy enough to compute the minimum 
degree of superheat required to nucleate a bubble from the 
surface. As can be seen from Fig. 3.5, as the cavity size 
increases, required minimum degree of superheat decreases. 

This evidently means that a given degree of superheat should 
activate a particular size of the cavity and theoretically 
speaking all other larger cavities existing in the surface. 

This refers to a particular size range of cavities that must 
be activated for a given degree of superheat. Pool boiling 

experiments reveal that the cavities which are actually acti- . , 

o-'sc cArfe-v ft un #>e co.vuii.-d we- 

vated probably because of their inability zo entrap vapor. 

A 

Equation (3.18) can be solved for R for given values of 
(T W ~T ), 6^, © and boiling liquid properties. This solution 
will give only two positive real roots which will comprise 
the size range of active sites. The size range computed using 
the proposed TD in the nucleation model is given in Table 3.3, 
along with those due to Hsu (1962) and Han and Griffith (1965). 

The cavities near the lower end of the size range nucleate 
while those near the upper and do not. Although the predicted 
range covers a wider spectrum of cavities, it is still pretty 
hard to say whether or not these large cavities will nucleate 
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or serve as stable nucleation sites. Although the experimental 

evidence of the activity of smaller sites is available , yet the 

same questions could, be raised about- them. Possible explanation 

might be that the smaller cavities may entrap vapor more easily 

and may generally have the geometry (L/D ratio etc.), to act as 

c 

stable sites. In this context, nucleation of cavities in the 
upper size range may be decided by other factors as well and not 
the degree of superheat and thermal layer thickness alone. 

Solution of equation (3.19) for given R c , (T -T^), 6^. and 
0^, results in a single cavity size which must nucleate. Any 
cavities larger than this will nucleate if other requirements 
(geometry etc.) are satisfied. It may be noted here that this 
equation can be used only for subcooled pool boiling. 
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DYNAMIC MODEL FOR BOILING NU CL BAT I ON STABILITY 

4.1 INTRODUCTION 

During boiling, vapor bubbles nucleate and grow from 
prescribed sites in the surface, called as active sites. Those 
sites which cannot trap vapor within them are called inactive 
sites. The steps in the life cycle of a bubble are shown in 
Fig. 4.1 which are based upon the experimental observations of 
Kant (1983) on single nucleation sites'. His photographic data 
showed the movement of liquid-vapor interface within the cavity 
and the presence of some liquid at the bottom of the cavity. 

As soon as a fully grown bubble lifts off from a cavity, 
the cold liquid partly from the thermal layer and the rest 
from the bulk rushes into the cavity and penetrates upto a 
certain depth. The temperature of this liquid, T , may be such 
that T, < T^ < T . The initial penetration is assumed to take 
place in a negligibly small amount of time. At the same time 
it is assumed that a very small amount of liquid left at the 
bottom of the cavity, at the vapor temperature when a bubble 
departs, cools down almost to the saturation temperature of the 
liquid due to flow of bulk liquid into the cavity at the moment 
of departure. Thus an elongated vapor bubble is entrapped 




Fig. 4.1 Diagram of proposed nucleation cycle showing the 
upper and lower l-v interface positions. 
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inside the cavity between the liquid at the top and that at 
the bottom. Owing to the curvature of the 1-v interfaces? the 
vapor trapped inside the cavity is at a higher pressure than 
that of the liquid. Since the temperature of the liquid in 
the Upper Liquid Control Volume (ULCV) , T , and that in the 
Lower Liquid Control Volume (LLCV) , T , are less than the 
saturation temperature of this vapor, condensation begins to 
occur. The upper 1-v interface recedes into the cavity and 
the lower 1-v interface rises up inside the cavity as the 
condensation continues. At the same time, since the surface 
temperature, which has to be higher than the vapor temperature 
for the boiling to occur on a metallic surface, conducts heat 
from the cavity walls into the ULCV and the LLCV, which 
increases the liquid temperatures. At the same time the 
condensation heat transfer across the two 1-v interfaces, 
further increases the liquid temperature inside the control 
volumes. The condensation phase ends when the temperature in 
both the control volumes reach the vapor temperature, T . The 
time elapsed during this phase is the condensation time, t . 
The wall temperature during the process is assumed to be at 
a constant value, T . 

Once the liquid temperatures have reached the vapor tem- 
perature, any additional heat transfer from the cavity walls 
goes to evaporate the liquid at the 1-v interfaces. This 
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introduces the evaporation phase which has been split into two 
phases, namely, first evaporation phase and second evaporation 
phase . 

During the first evaporation phase, the heat supplied from 
the cavity walls to the liquid in the ULCV and LLCV starts the 
evaporation of the liquids at the two interfaces- During the 
evaporation, the upper 1-v interface starts to rise toward the 
cavity mouth and the lower 1-v interface starts receding 
towards the bottom of the cavity. The first evaporation comes 
to an end when the upper 1-v interface just flushes the surface. 
By this time, it is assumed that the lower 1-v interface is at 
its lowermost position. The time elapsed during this phase is, 
t , as shown in Fig. 4.1. 

During the second evaporation phase the upper 1-v interface 

just protrudes out of the cavity and appears as a nucleus V' ■_ 

at the mouth of the cavity. During this period, the lower 1-v 

interface remains at the bottom of the cavity. The time elapsed 

during -this phase is t . 

e 2 

The sum (t +t ) is the total evaporation time. The sum 
e l e 2 

of the condensation and evaporation times (t +t ) is the waiting 

period, t , for a nucleating bubble. This bubble nucleus grows 
w 

at the site to its maximum diameter during the growth period, t . 

y 

The sum (t +t ) is the period of the bubble cycle which repeats 

W Cf 


itself . 
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If the liquid temperature never reaches the vapor tempera- 
ture? all the vapor will continue to condense, leading to a 
complete collapse of the cavity. 

The various steps shown in Fig. 4.1 are : 

a. Departure of the fully grown bubble from the surface and 
initial penetration of the liquid into the cavity. 

b. Condensation at the liquid-vapor interfaces. 

c. End of condensation. 

d. Evaporation at the liquid-vapor interfaces. 

e. End of first evaporation phase. 

f. End of second evaporation phase. 

g. Growth of the bubble at the cavity mouth to its maximum 
diameter. 

Based upon the above considerations, the dynamic model of 
Kant (1983) and Kant and Weber (1986a, 1986b, 1986c) for the 
nucleation of a bubble in HPS of a liquid is extended here to 
predict the waiting time, t . 

4 ,2 ASSUMPTIONS MADE. IN THE MODEL 

The nucleation cavity is cylindrical in shape. 

The 1— v interfaces retain their spherical shape during 
motion inside the cavity. 

The pressure just outside the cavity, p^ , is equal to the 

o 

system pressure. 
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4. At zero time, there is an initial penetration of the 

liquid into the cavity and a small amount of liquid is 

present at the bottom of the cavity at the saturation 

temperature , T , 
s 

5. The pressure and density of the trapped vapor within the 
cavity remain constant during the nucleation process. 

6. The vapor temperature, T , remains constant during the 
bubble nucleation. 

7. The temperature of the surface and that of the cavity 

walls, T , is a constant, 
w 

8. There are no thermal gradients existing within the vapor 
inside the cavity and hence the temperature of the inter- 
face is the same as the temperature of the vapor. 

4.3 CONDENSATION PHASE 

The initial penetration is assumed to take place in a 
negligibly small time. Condensation starts immediately at the 
1— v interfaces. 

4.3.1 Mass Balances 

During the condensation phase, as the upper 1-v interface 

( ULVI ) moves info the cavity some liquid from the bulk flows 

into it. The lower 1-v interface (LLVI) rises up inside the 

cavity as the condensation proceeds and some vapor is condensed 

at the interfaces. Let rL , m and A represent the mass flow 

*1 C 1 c 2 
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rate of liquid flowing into the cavity and the rates of conden- 
sation of vapor at the upper and lower interfaces respectively. 

Let IVU , iVu and M represent the mass of liquids and vapor 
1 *2 V 

respectively at any time t within the ULCV, LLCV and vapor 


control volume (VCV) shown in Fig. 4.2. 

The mass balances over the liquid and vapor control 
volumes yield. 


• ♦ 

rru + m 
n i C 1 



(4.1) 


dM 

-<% + V " ^ 

'■2 

% = ar - 


(4.2) 


(4.3) 


hU , itU and M are given by 
12 


tcd: 


z i~ V S c ] P 1 


(4.4) 


it d: 


NL, = [-5* (L-z 2 ) - V sc ] P 1 


(4.5) 


M v = [~4~ ( z 2~ z l^ + 2V sc^ 


V 


(4.6) 


where, V » the 

Sv# 

and Chilton, 1973) 


volume of the spherical cap is given by (Perry 




Fig. 4.2 Liquid and vapor control volumes 
within a cavity during the conden 
sation phase . 
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(1-sin©; (2+sinG) 


(4.7) 


24cos 6 


During the movement of the two 1-v interfaces within the cavity, 
the contact angle hysteresis may be disregarded and hence V 

s c 

may be assumed to be constant. The densities, P v and P^ , of 
vapor and liquid, respectively, are taken to be constant at the 
boiling point of the liquid. 

Adding equations (4.1), (4.2) and (4.3) and using equations 
(4.4), (4.5) and (4.6), yields 


dz, dz 2 

T~ ( p l “ p v } ^dt" ~ dt"^ 


(4.8) 


Adding equations (4.1) and (4.2) and using equations (4.4) and 
(4.6) yields 


nt, - = 


7 xQ dz. dz 0 dz. 

-4 s 1 5 sr + p v ( ar - ht^ 


(4.9) 


The condensation heat transfer across the LLVI is given by 


q ic 2 h ic ; ^ T v~ T i 


(4.10) 


where h. „ 
xc 


film heat transfer coefficient for condensation 
interfacial area given by 


(4.11) 


2(l+sin8) 
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If is the enthalpy of vaporization of the liquid, the mass 

flow rate of the condensate, m , into the LLCV is given by 


q. %\T h. (T -T„ ) 

rc 2 c ic v v 

C 2 ^fg 2( l+sin©) h- 


(4.12) 


Let be the average flow velocity inside the cavity in the 
ULCV, the mass flow rate, m „ can also be written as 


: , m , J 


-r 5 V 1 


(4.13) 


Substituting for m and m , from eauations (4.12) and (4.13) 

c 2 *1 

into equation (4.9) yields, 


dz, i 2h ic^ T v" T 1 J dz 2 

1 = ~ — — * [ P 1 V, — — - - P v dF"] 

dt (P 1 - P y ) (l+sin©) h fg 


(4.14) 


Again from equations (4.2) and (4.3), we get 


= dF + dF 


(4.15) 


considering the condensation heat transfer across the ULVI , the 
mass flow rate of the condensate can be written as : 


qic i 1t ° chic(Tv ~ Tl 

: 1 h fg 2 (l+sin©) h f 


(4.16) 


Substituting for M 1 , M y and m from equations (4.5), (4.6) 
and (4.16) into equation (4.15) yields, 
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dZ r 

dt" 


p dz, 

'1 ■" v' Ti+sInQ) h^g v dt 


(4.17) 


4.3.2 Momentum Balances 


The momentum balance over the ULCV shown in Fig. 4.2 


yields 


ntf; ittf: 

(P 1T V ~ + ( “4- 2 1 - V sc> »1 9 


2 

tcD dz-, 

+ 3t~ z l> V 1 - 1lD c 2 l T w 


r, UD?Z, 

2™ [ — £__i _ u 1 p v 
dt L 4 V sc J V 1 


(4.18) 


or 


dV-j 

dt” 


= g + 


[*i - 


4V 


C(p„ “P« ) - — jjj— (4.19) 


sc 


TCU 


K 


ii «i 2 


2 -4 


where 


p n 


liquid pressure just within the mouth of the 
cavi ty 


p . = liquid pressure just above the ULVI 

1 2 


t = wall shear stress exerted by the wall on to the 
w 

liquid within the ULCV • 
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a. Pressure loss due to sudden contraction 

There is a pressure or head loss due to sudden contraction 
involved when the liquid flows into the cavity from outside. 

It is given by 

\p- 

h^ c = 2g me ^ er °f fluid (4.20) 


where C c is the coefficient of contraction which can be taken to 
be approximately 0.5. This head loss can be expressed in terms 
of a pressure loss as follows : 



where C-^ 


= 1.333224 x 10“ 


c r 

<2i )( 



(4.21) 

(4.22) 



= density of mercury at room temperature. 


Therefore, the pressure p . can be written as 

1 1 

P„ = p - C,V? (4.23) 

'll 1 1 1 

o 


b. Determination of p, 

12 


The pressure, 
equation : 



can be related to p y by the Laplace 



(4.24) 
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where a is the surface tension of the liquid at the boiling 

point and the radius of curvature of the spherical interface 

which is given by cos 9 = R . The pressure, P. , can then 

n 2 


be written as : 



4a c os 6 
D 


c 


(4.25) 


Substituting the values of p and p from equations (4.23) and 

n l n 2 

(4.25) into equation (4.19). 


dVj 

dt" 


= 9 + 


4V 

0 


[( Pl _ - P v - 


, 4a cos 0\ ^1 Wt 
+ n > ~ n J. 


tcD 


'I 


(4.26) 


Determination of wall shear stress 


The velocity distribution inside the cavity is taken to be 
parabolic as in the case of Hagen-Poiseuille flow. The velocity 
v in the z-direction at any radius r is given by 

v = 2V 1 [1 - (f-) 2 ] (4.27) 

l c 

From Newton’ s law of viscosity 

dv 2 

T w ~ ^dr ^r=R c 

8V x p 

T w = ~ 


(4.28) 
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Substituting into equation (4.26), we get 


dV-, 

dt" 


g + 


[z 


4V 


1 - 

X %D Z 1 


f f rv — D p tP- . COS S'! 

UP 1 P v C 1 V 1 + ~D } 


32 


V z i. 


D" 


(4.29) 


Consider now the LLCV. Unlike the ULCV, there is no liquid 
inflow into it. The 1-v interface will tend to rise up in the 
LLCV with a velocity of V 2 (-ve velocity), only due to the con- 
densation at the LLVI . It is assumed here that the wall shear 
stress exerted by the wall onto the liquid acts only due to the 
movement of the 1-v interface. The momentum balance over the 
LLCV yields 

2 2 

~ ( p 1 4 " P 1 3 ^ T" + ^~"4~" ( l -z 2 )” V s J P 1 g + %D c T w ( ‘ L ~ z 2^ 

d %D c 

= at Hr ( L “Z 2 )-V S c]P 1 V 2 (4.30) 


or 


dV r 

dT 


g - 


[ (L-z 0 )- 


4V 


sc 


it D~ 


] P> 


[(P, -P. )- 

» A 1 Q 


4(L-z 2 ) t 
D 


w . p v 

P 1 v 2 dt J 
(4.31) 


The wall shear stress is given by 

8V 2 p 

T w = -"T" 
c 


(4.32) 
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Substituting equation (4.32) into equation (4.31) one obtains 


dV 9 

dt~ = ^ 


4V 


sc- 


[(L-z 9 )' 
32V 0 p(L-z 0 ) 


C(P -P ) 

14 13 


+ 


D 4 


dz 0 

- P-,v 2 ^-] 


(4.33) 


Determination of p and p 

>3 M 

The pressure, p^ , is assumed eaual to p., , obtained from 

3 ‘ *2 

equation (4.25). 

The pressure, p. , is given by 

1 d 


= p^ + [z x + (L-z 2 )] g (4.34) 


4.3.3 Energy Balance 

Considering the ULCV shown in Fig. 4.3. 

Sensible heat of the liquid flowing into the cavity 

= m c ( T -T.) (4.35) 

"1 x 1 o b 

Sensible heat of the condensate going into the ULCV 


“c, WV 


(4.36) 


Heat transfer from the vapor to the ULCV through the interface 

(ULVI > = A, hi (VT, ) 

^ X 


(4.37) 



— Hfg 


m C 2 hfg 


tics for the liquid control 
g the condensation phase. 
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Transfer of heat from the cavity walls 
~ ^wli^ Nvli^w~^1 ^ 

The expansion work due to expansion of the ULCV 

—.2 

l a 

= A n f , £ _ 1 

J P 1 2 dt L 4 z l^ 

Change in internal energy of the ULCV 

= 3t u l> = 3 t£\ 

The heat balance over the ULCV yields, 


(4.38) 


(4.39) 


(4.40) 


% c 1< W + '"c 1 C l( T 1 1 - T b) + Aic h ic<V T 1 1 ) 

„2 

i j % D 

+ 'W h wli (VTl l ) = 5' \ ™ [ ~*" ^ 

+ C 1 < T 1 -T b >3 ( 4 - 4 D 


' utting Ajl 4 C (T V -T n ) + (T^-Tj = A i h ic (T v -T 1 ^ 


A , . = 7i D z, 

wlij_ c 1 


(4.42) 


(4.43) 


Using equation (4.42) and substituting the values of A^, 
A w n , and from equations (4.11), (4.43), (4.13) and 
(4.4) into equation (4.41) and simplifying yields 




Considering now the LLCV as shown in Fig. 4.3. 
Sensible heat of the condensate going into the LLCV 

= Vi (T 1„-V 


(4.45) 


Heat transfer from the 
( LLVI ) 


vapor to the LLCV through the interface 


= \ h ic ( V t i 2 ) 

Transfer of heat from the cavity walls, 

A wli 2 h wli^ T w~ T 1 2 ^ 

The expansion work due to the expansion of the LLCV 

n 2 

1 d 

= J p l 3 Ht (l-z 2 )] 

Change in internal energy of the LLCV 


(4.46) 


(4.47) 


(4.48) 


[m i 2 c i <V T * )] 


(4.49) 
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The heat balance over the LLCV yields . 


“c 2 c 1 (T V T b) + A i 4c ( T v- T 1 2 ) + Awli 2 h wli4 w -T l2 ) 
= dt ^ M 1 2 °1 ^ T 1 2 ~ T b^ + J P 1 3 3t ^~4~ ^ L ~ z 2 ^ 


(4.50) 


Putting m c1 <4 -T fa ) + 4 h ic<V T 1 > 


= A. h. (T -T. ) 
l ic v 


(4.51) 


and 


\li 2 ~ T5 D c^ L ” z 2^ 


(4.52) 


Using equation (4.51) and substituting the values of A^, \jij 
and from equation (4.11), (4.52) and (4.5) into equation 
(4.50) , one obtains 


dt 


^ic ^ T v 


Tl 0 > 


4V 


(1+sin 9) P. c, [(L-z 9 ) - — 
11 tiD, 


4h wli (L - Z 2) ( V T 1 2 ) 


4V 


D p c [(L-z 2 ) —I 

c > * z %YT 


+ 


4V 

[(L-z 2 ) f 2 ] 

_ c 


1 3 dz^ 
[(T. -T h ) + ] 3T 


12 


Jp., C A 
1 1 


(4.53) 
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4.3.4 Initial Conditions 


iho liquid from the bulk rushes towards the site, immedia- 
tely o£jNr a bubble departs from that site. This liquid pene- 
trates to a certain extent into the cavity in a negligibly small 

time. At the moment of bubble departure i.e., at t = 0, there 

* 

is a small amount of liquid present at the bottom of the cavity. 
Both the ULVI and LLVI are spherical cap shaped whose height is 
given by (see Fig. 4.2) ; 

S = 2 D c (sec © - tan 9) (4.54) 

The liquid that penetrates into the cavity is assumed to 
have mean velocity of V^. Here a value of 1.0 m s” 1 is taken 
in conformity with the observations of Anton and Vekas (1973). 

The velocity of the LLVI and hence the induced velocity of the 
liquid in the LLCV, V ^ s is assumed to be zero at t = O. 

The liquid that rushes into the cavity and enters the ULCV 

may come partly from the bulk and partly from the thermal layer. 
The temperature of the liquid within the ULCV, may be such 

that at t = 0, Tfo < T 0 < T w . 

^The presence of this liquid can be clearly seen from the photo- 
graphic data of Kant (1933). It may either be present due to 
the flow of bulk liquid to the cavity bottom at the moment of 
departure or it may be present to maintain a contact angle 9 
for the given surface-liquid combination. 
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The liquid at a temperature, T , present in the LLCV is 

2 

assumed to be cooled to T g , possibly due to the penetration of 
the bulk liquid into it at the moment of bubble departure (i.o., 
at t = 0) . 

The positive direction is taken into the cavity and the 
coordinates and measure the distances of the ULVI and 
LLVI at any time t from the mouth of the cavity. 


The initial conditions for the condensation phase can be 
summarized as follows : 


At t = 0 


Z 1 = 


V, = 


Vo = 


'll 

‘12 


L-S 

\ 

0 

T 


= T 


0 < V Q < 1 


T < T < T 
b — o — w 


(4.55) 


Condensation phase is considered to be complete as soon as 
Tj = T^ = T y . If any of the liquid temperatures reach the 
value of T y earlier, it is taken to bo constant until the other 
liquid temperature reaches T . Equations (4.14), (4.17), (4.29), 
(4.33), (4.44) and (4.53) were solved numerically subjected to 
the initial conditions given by equation (4.55) to determine the 

c* 


condensation time t 
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4.4 FIRST EVAPORATION PHASE 


When the condensation phase is complete, any additional 
heat transfer from the cavity walls will mostly go to evaporate 
the liquid at the ULVI and LLVI causing reversal of their motion. 
The first evaporation is assumed to be over when the top of 
the ULVI is even with the surface, i.e., z ± = S. By this time, 
the LLVI reaches its maximum value, = L-S. The liquid 
present in both the ULCV and LLCV are at the vapor temperature 

T . 

v 

4.4.1 Mass Balances 


During the upward motion of the ULVI, some liquid will 
flow out of the cavity. During this period, the LLVI will move 
downwards until it almost reaches the bottom of the cavity. The 


- m - rii 
12 


m 


m 


+ m 


and LLCV shown in Fig . 

4.4, yield 

dM . 


h 

dt 

(4.56) 

dM 


V 

" dt 

(4.57) 

dM 


1 ? 
dt 

(4.58) 


where m and m denote the rate of liquid evaporation at the 
e l ®2 

interfaces and m. rate mass outflow from the cavity. Adding 

‘2 

equations (4.56), (4.57) and (4.58) and using equations (4.4), 
(4.5) and (4.6) yields 




Fig. 4.4 Liquid and vapor control volumes within 

a cavity during the first evaporation phase 
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-m. = 


U~ pj [ 


dz^ dz 2 . 


*1 2 4 V L dV 


(4.59) 


4.4.2 Momentum Balances 

The momentum balance over the upper liquid control volume 

shown in Fig. 4.3 yields 

2 2 
tcD'T tD 

-(P i -P, ) + ("^ z - V )p„ g + itD z. t 

<2 *]_ 4 4 i sc' -] ^ c 1 w 


7rD c dz l ri 

= ~ 4 ~ P 1 <- V l> dT + I 


[ C _ v J p v (4.60) 


dt T dt L 4 A l " v S c J ^ V 1 


Rewriting the above equation 


4V ^ P 1,“ p 1 2 ' + "D ^ 
r _ sc 1n '1 l 2 c 

tZl ’7 dT' j i 

C 


(4.61) 


a» Pressure Loss due to sudden expansion 

There is a head or pressure loss h due to sudden expan- 

le 

sion involved when the liquid flows out of the cavity mouth, 
and it is given by 


C o~“ meter of fluid 
e 2g 


(4.62) 


where C e is the coefficient of expansion approximately equal to 
unity. This can be written as 


c 2 v^ 


(4.63) 
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whore C 2 = 1.333224 x 10 5 (^g)(£l~) 

’ P Hg 

The pressure, , can be written as 

Therefore, p = p + 0 o V? 

n i 0 z 1 


(4.64) 


(4.65) 


b. Wall Shear Stress 


The wall shear stress for the Hagen-Poiseuille flow using 
the velocity distribution given by equation (4.27) is given by 


^ T w^HP 


8V. p 


(4.66) 


Substituting the values of p. , and t from equations (4.65) and 

I i w 

(4.66) respectively, one obtains 


dV. . 9 32V. pz, 

= g + i_. [(p +C 2 W-p ) - — 4-±] (4.67) 

Cz. fS] p 0 2 D c 

c 

The momentum balance over the LLCV shown in Fig. 4.3 yields, 

_2 2 
it LT uDt; 

~( p i o “ P 1 A “ 4 ™ + f“ 4 “ (L” z 2 ^" V sc-^ P 1 9 " ^ d c ( l - z 2 ^ T w 

O *T 


, rcDt: 

dt ^T* ^ L_z 2^“ V sc^ P 1 V 2 


(4.68) 


Substituting t = 
gives 




in the above equation and simplifying 


32V9M-(L-Z9) 

L -4t r [( W + — 

[ ( l -z 2 ) - KD 2 ]P n 
c 


dZ9 

- 6 V 2 

(4.69) 
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4.4.3 Energy Balance 

During the evaporation phase, the liquid temperatures, 

T , and, T , in the ULCV and LLCV are already at the vaoor 
1 1 1 2 

temperature, T^, and there is no heat transfer across the 1-v 
interfaces. For computing the energy quantities for the control 
volumes, is taken as the reference temperature. 

Consider the ULCV as shown in Fig. 4.5. 

Sensible heat of the liquid leaving the ULCV 

■ % p ( W ' 4 - 70) 

Heat of vaporization of the mass leaving the ULCV 



Transfer of heat from the cavity walls 


(4.71) 





iquid control 
aporation phase. 




63 


Change in internal energy of the ULCV 


dM 


AUl= ft u) 


= M 


du 


'll 


+ u 


dt 


dT 

WL c. xr + c -i ( T„“T k ) 


dM, 


“»1 ^ dt 


'1 v v b' dt 


(4.74) 


=0 as T v = constant. 

Energy balance over the ULCV yields 

- m 1 c 1 (T v -T b ) “ m e;L h fg + ‘Vl^ h wli (T w“V 

x 71 dz 1 

= C 1 ( T v -T b ) » d % + J P 1 2 T" dT 


where 


dt 


A , . = TtQ.Zi 

rt wli b c 1 


(4.75) 


(4.76) 


and using e 


quations (4.57) and (4.53) m can be written as 


,,, dM, 

dM v • _ J2 

= ‘dPT “ m en “ dt + dt 


i ~~ -2 

equations (4.5) and (4.6) gives 


which upon using 

;D S 

T” L ^v dt 


nrf d z x . v dz 2 1 

-m o = + — L p v dT + (p 1 ' P v } d ^ 
"1 


(4.77) 


and M, from equa- 


Substituting the values of m i 2 ’ 

tions (4.76), (4.59), (4.77) and (4.4 ) respectively, into 
equation (4.77) and simplifying yields 
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* ^ h fg^ P 1 " P v> 


+ P v c 1 ^ T v~ T b^ ~ p v h fg^ 


dz n 4z, h , . 

- c,(p - P V )(T V -T b )] ( W } (4.73) 


Consider now the LLCV. 

Heat of vaporization of the mass leaving the lower control 


volume 


nr h .p~ 

e 2 fg 


(4.79) 


Heat transfer from the cavity walls to the liquid 


\/li 2 ^wli (T w T v> 


(4.80) 


The flow work due to the contraction of the LLCV 


1 JWJ J 

1 n --A SL. ( T-Z ) 

J P 1 3 4 dt U Z 2> 


(4.81) 


Change in internal energy of the lower control volume 


dt 


(M„ u) 


dT v fT T x 
f S 2 C 1 dt + c 1^v Hr dt 


C 1 O v -T b ) gt- 


(4.82) 
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Energy balance over the LLCV yields 


. 1 tcD , 

" m e 9 h fg + \jli~ h wli ^ T w“ T v^ = J p -i “ 4 “ dt ^ L " z 2^ 


diiL 


+ c , (T -Tj 


1 ' v b' dt 


(4.83) 


Substituting for m e , Mj and A wlj _ from equations (4.58), (4.5) 

2 2 2 

and (4.52) into equation and simplifying yields 


dZr 

d'F 


r - 

[ p v h fg -T- p i c i (V T b>l 


4 W L - z 2>< T w -T v ) 

[ ■] 


(4.84) 


4.4.4 Initial Conditions 

The conditions at the end of the condensation phase are the 
initial conditions for the first evaporation phase. 


At t = t 


Z 1 “ z cl 


z 2 " z c2 


V, = V 
1 me. 


(4.85) 


V 2 ^ mc 2 

Equations (4.78), (4.84), (4.67) and (4.69) were solved, numeri- 
cally subject to the initial conditions given by equation (4*85) 

to determine the time t corresponding to =S. 

e l 
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4.5 SECOND EVAPORATION PHASE 

During this phase, it is assumed that heat is transferred 
in the axial direction as well as radially from the cavity walls. 
For the vapor nucleus to appear above the surface, an addi- 
tional time t must elapse before becomes equal to zero. The 
liquid present at the bottom of the cavity is assumed to play no 
role during this phase. 

The surface area of the spherical segment above the heating 
surface, as shown in Fig. 4.6, -is given by 


A im = 2ltR i (S - z l> 


(4.86) 


Assuming that the liquid temperature at the cavity mouth 
at t = t equals the wall temperature T , the rate of evapora- 

Q w 

tion at the ULVI is 


7rD c z l h wli^T'w“I v ) + A^h^ (T w -T v ) 


(4.87) 


where 


h. = film heat transfer coefficient for the evaporating ULVI 
le 


(assumed equal to h- c ) 


.0 

Substituting for A^ m from equation (4.86) and writing = cos q 
into equation (4.87), yields ;| 


AlVU [(h 


h. Sh- _ 

_ 10 ■ h - JL. — I 

wli cos Q 1 cos© 


(4.88) 





Fig. 4.6 l - v interface protruding from a cavity 
during the second evaporation phase. 
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Using equation (4.77) and assuming (dz 9 /dt) = 0, yields 


dz 

Ht 


1 _ 


4 (T w~ T v) j- / , * Sh ie 1 

D c v h fg h wli ~ cos Z 1 + coT © J 


(4.89) 


The momentum equation for the second evaporation phase is the 
same as equation (4^67) . 

Initial Conditions 

The conditions at the end of the first evaporation phase 
are the initial conditions for the second evaporation phase. 



4.6 PREDICTIONS 



(4.90) 


The model equations for each of the three phases, developed 
in the preceding sections were solved numerically to predict 
the waiting time of a bubble for given boiling conditions. The 
method of solution of the equations is described in Appendix C. 

The equations were used for both subcooled and saturated 
pool boiling. ■ 

4.6.1 Application of the Model 

The solution of the equations involved the values of h ie and 
h if which were calculated using equations (D.l) and (D.2) . The 
multipliers and iV^ were chosen for water to agree with the 
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available frequency data of Shoukri and Judd (1973). The 
waiting time wa^ taken as lb/ of the bubble cycle period. In 
accoroance wi'&h the measurements of Shoukri, the contact angle 
for water boiling on a copper surface was taken to be 22°. For 
isopropanol, the contact angle was taken as 10° and Mg was 
related to that for water through a ratio of thermal conducti- 
vities (Kant, 1983) as 

M 2 = (k, /kip 1 - 25 M 2r (4.91) 

where Mg and k-, represent the multiplier and the thermal con- 
ductivity for isopropanol and M 2r and ki r are the same for the 
reference liquid which was taken as water here. The values of 
and M 2 for both the liquids are given in Table 4.1. 

Table 4.1 

Multipliers for h^ c and h ^ 


Boiling 

Liquid 

M x 

Mg 

Water 

IxlO -4 

5xl0“ 3 

I sopropanol 

-4 

1x10 

-4 

7x10 
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4.6.2 HPB of Water and Isopropanol 

The model equations were solved numerically for both 
water and isopropanol. The various parameters were computed 
as a function of time. The properties used for the various 
liquids are tabulated in Appendix A. The solution of equa- 
tions predicted the waiting time as a function of cavity 

diameter , cavity-depth, degree of superheat. (T -T ) . and the 

w s * 

degree of subcooling, (T -T^). 

The predictions for two different cases of subcooling 
and different degrees of superheat for water are given in 
Tables 4.2 and 4.3. Table 4.4 gives the predictions for the 
saturated boiling of isopropanol. These tables also give the 
computed values of S, p v , T v , h ic , h wli , z cl and z c2 - 

The positions of the ULVI and LLVI at different times 
are shown in Fig. 4.7. It is seen that the ULVI penetrates 
into the cavity to its maximum value, = z c j_» within a 

fraction of a millisecond. The reversal of velocity from 
positive to negative takes place as soon as = z is 
reached and then slightly decreases. The ULVI 

stays at this position, more or less, throughout the condensa- 
tion phase, i.e., until the liquid temperature T ^ reaches the 
vapor temperature, T . Once this condition is reached the first 
evaporation phase begins during which the ULVI rises up In the 
cavity to a value z^ = S • This phase lasts for a fraction of a 
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millisecond. The second evaporation phase begins as soon as 
z becomes equal to S and the ULVI rises further up towards the 
cavity moAth until z-^ becomes equal to zero and the vapor 
nucleus protrudes out of the cavity. 

The LLVI during the whole cycle is more or less at the 
bottom of the cavity where = L-S. 

During the condensation phase, the velocity, V^, drops to 
a small value in a fraction of a millisecond and turns negative 
and shows a slight' increase in the value. The velocity, V 2 , 
rises upto a few milliseconds in the negative direction and then 
stays constant at a negligibly small value. 

The variation of temperatures Ti and Ti are shown in 

1 2 

Fig. 4.3. The temperatures of the liquid present in the ULCV 

# and LLCV, increases steadily during the condensation phase 

until they become equal to the vapor temperature. The liquid 

in the ULCV takes a longer time to attain the vapor temperature 

than the liquid in the LLCV. The temperatures T -j and T-j 

1 2 

remain at T during the evaporation phases. 

4.6.3 Effect of Model Parameters 
a. Cavity Diameter 

The predicted waiting time increases as the cavity dia- 
meter, D increases. This is because more liquid from the bulk, 
c 

flows into the cavity which then takes a longer time to heat up 
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to the vapor temperature. Also as D increases; p , T , h. 

0 V V 3* 0 

and h ^ decrease which result in reduced rate of heat 
transfer into the cavity 9 and larger depth of penetration 
thereby increasing the waiting time. 


As shown in Tabid 4.2, theoretically predicted value of 
the waiting time is in good agreement with the experimental 
value . 


b. Degree of Superheat 

Tables 4.2, 4.3 and 4.4 also show the computations for 
various values of degree of superheat. An increase in its value 
increases the heat transfer rate and thereby results in a lower 
waiting time. It can also be seen that the depth of penetra- 
tion is not significantly dependent upon the degree of super- 
heat. 


c. Degree of Subcooling 

! 

A larger degree of subcooling increases the predicted 
waiting time, but the depth of penetration is not affected by ; 
it (see Tables 4.2 and 4.3). | 

4.6.4 Stability of Nucleation Sites ! 

The model equations presented earlier relate all the rele-: 
vant parameters of boiling nucleation phenomena. Computations | 
for given boiling conditions as shown in Tables 4*2, 4*3 and 4 - 
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•tv size which controls Zj. 

show clearly that it is only the ^ nucle ation site, 

' t h pthor or not a cavity i= , v < V . 

and z 2 - -hschci ^ ^ . if (l/D c ) < * 7 v - * 

* U1 depend UP0 " ^ V " fined n/with the boiling 

and the cavity 9= 




CHAPTER 5 


CONCLUSIONS 

5.1 MODELING OF THERMAL LAYER THICKNESS 

All the relevant parameters of N.B. were considered in 
developing an expression for The dimensionless groups 

obtained were combined in an arbitrary manner to reduce the 
number of data points required to evaluate the indices 
(i = 1, 4). The results obtained are in good agreement with 
the values of 6 ^ reported in literature. 

Unfortunately, there is no reliable experimental data on 
6 ^ because of which only estimated values of 5 ^. for various 
liquids were used. The apparent limitation of the proposed 
expression is that some experimental data on 6 ^ is definitely 
required to evaluate the indices . But the dimensionless 
groups obtained could be combined into one or 2 groups if a 
very limited data are available. 

5.2 BUBBLE NUCLEATION SUPERHEATS 

Kant's (1983) bubble nucleation model was extended to 
incorporate a proposed TD given by equation (3.14). Substitu- 
tion of various parameters in equation (3.18) and the constants 
given in Table 3.2 resulted in Nucleation Superheats somewhat 
higher than those predicted by Kant, and naturally in better 
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surface and evolve criteria for the stability of these sites. 

Kant's experiments clearly exhibited the presence of some 
liquid at the bottom of the cavity and thus in an active 
nucleation site, there was liquid present at the top and the 
bottom and a vapor bubble was stuck in the middle. This pro- 
vided the l~v interfaces at the top and bottom (ULVI and LLVI). 
The model proposed here for the first time took into account 
the condensation and evaporation at both of these interfaces 
to predict the waiting times of bubbles which were found to be 
in agreement with the available experimental data of Shoukri 
and Judd (1978) . The predictions of the model also showed that 
the diameter and the depth of a cavity basically determined the 
stability of a nucleation site. Other variables like degrees of 
superheat and subcooling were not found to be important. 

The model was applied to subcooled and saturated pool 


boiling . 
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APPENDIX A 

PROPERTY DATA 

The property data at the Boiling points were calculated 
as follows : 

1. The density of a pure liquid, , at the B.P. was 
calculated (within + 3 percent) using Benson’ s equation (Perry 
and Chilton, 1973) 

P 1 . = P (1.981 + 0.422 log pj (A.l) 

where 

-3 

P = liquid density at the critical point in g-mol cm 
p = critical pressure of the liquid in atmospheres 

V 

absolute 

The liquid density was also calculated using Tyn and Calus 
method (Reid et al. 1977). 

The method is used for estimating liquid molal volumes, 

V b , a t B .P . 

V, = 0.285 ’g 1 * 048 (A. 2) 

3 ~i 

where V = critical volume in cm g-mol taken from the 
c 

Property Data Bank of Reid et al. (1977). 
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2. The density of the vapor, P , was calculated using the 
Ideal gas equation at the normal B*P» and atmospheric pressure. 




(A. 3) 


where p = atmospheric pressure, 1 atm 
T = temperature, K 

R = universal gas constant = 82.04 J kg-rnol ^ K ^ 


3. The dynamic viscosity of a pure liquid, at the normal 
8.P. was determined using the equation suggested by Van Velzen, 
Cardozo and Langenkamp' s method (Reid et al., 1977). 


log = VISB(i - (A.4) 

where , 

VISB , VISTO = constants related to structure obtainable 
from the Data Bank of Reid pt al. 

T = temperature, K 
\x ^ = viscosity, centipoise. 


4. The surface tension, a, at the normal B.P. was estimated 
within « 5 percent using Walden’s rule (Perry and Chilton, 


1973) 


cr =s 


h fQ P 1 


364 


(A. 5) 
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where a = surface tension, dynes cm -l 

h fg = entholpy of v apori2ation ? cal g _ mol ~l 
P 1 = liquid density, g crn ~3 

5. The thermal conductivity of a pure liquid, ^ , at the 
normal B.P. was established uslng Sato , s correlatlon glv6n ln 
(Reid et al., 1977). 

, _ 2 .64x10" 3 

N (A- 6) 

M ' ' 

where 

k 1 = thermal conductivity, cal cm -1 s -1 K -1 
M = molecul ar weight 

6. The liquid specific heat, c^ , of Pentane and Ether at their 
B.P. was estimated using Mis$ e n a rd ! s group contribution method 
(Reid et al. 1977). The actual values at the B.P. were cal- 
culated using linear interpolation from the available data. 


% 



Properties of some pure liquids at Normal B.P. 
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.APPENDIX B 


DIMENSIONS OF VARIABLES USED IN CHAPTER 3 

The dimensions in Ms,L,T,9 system are given below : 


[6 ] 

= L 

Cq] 

= mlt” 2 

[T -T ] 

u w S J 

= 9 

t o v ] 

= ML” 3 

CPU 

= ML” 3 

tV 

= l 2 t” 2 

[0] 

= mt” 2 

tv 

= L 

[ v 

= ML” 1 !” 1 

Ck 1 ] 

= MLT” 3 ©” 1 

[c 1 ] 

2 -2 -1 
= Li ^9 1 

[S R ] 

= L 




APPENDIX C 


SOLUTION OF FIRST ORDER ORDINARY DIFFERENTIAL EQUATIONS 

The system of ODE's developed in Chapter 4 were solved 
numerically using the NAG Fortran Routine D02BDF (X,XEND,N,Y» 

TOL , IRELA3 , FCN , STIFF , YNORM , W, IW ,M , OUTPUT , IFAIL ) . 

P ARAMETERS : 

X - initial value of the independent variable T 

XEND - final value of the independent variable 
N - number of equations 

Y - array of dimension N 

Y(l), Y(2), . Y(N) must contain the initial values 
TOL - tolerance for error control 

IRELAB - determines the type of error control 

= 1 for number of correct decimal places 
= 2 for number of correct significant digits 
= 0 mixed error test 

FCN - subroutine for evaluating functions. 

It should be provided by the user and be of the type 
SUBROUTINE FCN(T,Y,F) 

DIMENSION T,Y(N), F(N) 

FCN must be declared as EXTERNAL in the calling program. 
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STIFF 


YNORM 

W 

IW 

M 

OUTPUT 


I FAIL 


- For stiffness check 

=0 no stiffness check is made 

>0*0 stiffness value is computed at the end of each 
integration step 

< 0.0 only the first stiffness value is calculated 

- If YNORM - 0.0 then it has no effect 

>0.0 used as a bound on the computed solution 

- array of DIMENSION (N, IW) used partly for output 
and partly for workspace 

- second dimension of W, > 14 

- to call routine OUTPUT every M integration steps 

- subroutine used to print the values of its arguments; 
It should be supplied by the user and be of the type 
SUBROUTINE OUTPUT (X,Y,W, STIFF) 

OUTPUT must be declared as EXTERNAL in the calling 
program. 

- error indicator. 
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APPENDIX D 


CONVECTIVE HEAT TRANSFER COEFFICIENTS 


The adjusted values of the convective heat transfer 
coefficients h ic and h ^ were determined using the following 
equations 


h ic = hl Mjl 


(D-l) 


h wli= h w “2 

The interfacial coefficient, h^, was determined from an 
expression based upon the kinetic theory of gases (Kant, 1983, 
equation (A4.4)) as 



The coefficient, h , was taken as the heat transfer coefficient 

w 

for the laminar fully developed flow of a fluid through a tube, 
for constant wall temperature, given by 

h , = 3.658 (k,/Dj (D.4) 

W l u 

Other symbols have their usual meanings. 



